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Abstract — This work considers the problem of quickest detec- 
tion with N distributed sensors that receive continuous sequen- 
tial observations from the environment. These sensors employ 
cumulative sum (CUSUM) strategies and communicate to a 
central fusion center by one shot schemes. One shot schemes 
are schemes in which the sensors communicate with the fusion 
center only once, after which they must signal a detection. The 
communication is clearly asynchronous and the case is considered 
in which the fusion center employs a minimal strategy, which 
means that it declares an alarm when the first communication 
takes place. It is assumed that the observations received at the 
sensors are independent and that the time points at which the 
appearance of a signal can take place are different. It is shown 
that there is no loss of performance of one shot schemes as 
compared to the centralized case in an extended Lorden min- 
max sense, since the minimum of N CUSUMs is asymptotically 
optimal as the mean time between false alarms increases without 
bound. 
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I. Introduction 

The problem of decentralized sequential detection with data 
fusion dates back to the 1980s with the works of [1] and 
[2]. We are interested in the problem of quickest detection 
in an iV-sensor network in which the information available is 
distributed and decentralized, a problem introduced in [16]. 
We consider the situation in which the onset of a signal 
can occur at different times in the N sensors, that is the 
change points can be different for each of the N sensors. We 
assume that each sensor runs a cumulative sum (CUSUM) 
algorithm as suggested in [7], [11 ]— [ 14] and communicates 
with a central fusion center only when it is ready to signal 
an alarm. In other words, each sensor communicates with the 
central fusion center through a one shot scheme. We assume 
that the N sensors receive independent observations, which 
constitutes an assumption consistent with the fact that the N 
change points can be different. So far in the literature (see 
[7], [1 1]— [14]) it has been assumed that the change points 
are the same across sensors. In this paper we consider the 
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case in which the central fusion center employs a minimal 
strategy, that is, it reacts when the first communication from 
the sensors takes place. We demonstrate that, in the situation 
described above, at least asymptotically, there is no loss of 
information at the fusion center by employing the minimal one 
shot scheme. That is, we demonstrate that the minimum of N 
CUSUMs is asymptotically optimal in detecting the minimum 
of the N different change points, as the mean time between 
false alarms tends to oo, with respect to an appropriately 
extended Lorden criterion [5] that incorporates the possibility 
of N different change points. As an observation model we 
consider a continuous time Brownian motion model, which is 
a good approximation to reality for measurements taken at a 
high rate. Moreover, given a high rate of observations from 
any distribution, the central limit theorem asserts that sums 
of such observations are normally distributed and therefore 
the Brownian motion model is a plausible model for such 
situations. 

The communication structure considered in this paper is 
summarized in Figure Q] in which T, for i = 1 , . . . , TV 
denote stopping times associated with alarms at sensors Si 



i = 1, . . . , N, respectively. 

In the next section we formulate the problem and demon- 
strate asymptotic optimality (as the mean time between false 
alarms tends to oo), in an extended min-max Lorden sense, 
of the minimum of N CUSUM stopping times in the case of 
centralized detection. We then argue that this result suggests no 
loss in performance of the one shot minimal strategy employed 
by the fusion center in the case of decentralized detection. We 
finally discuss an extension of these results to the case of 
correlated sensors. 

II. The centralized problem 

We sequentially observe the processes {Q ;t > 0} for all 
i = 1, . . . , N with the following dynamics: 



(1) 



« 



dw 



(i) 



t < n 



[idt + dw t t > Ti, 



where /i > is known 0, {«4 } are independent standard 
Brownian motions, and the Ti's are unknown constants. 

An appropriate measurable space is 17 = C[0, oo) x 
C[0,oo) x ... x C[0, oo) and T = U t> oJ r t, where {T t } 
is the filtration of the observations with Tt — cr{s < 



i) 



£i )}. Notice that in the case of centralized de- 



tection the filtration consists of the totality of the observations 
that have been received up until the specific point in time t. 

On this space, we have the following family of probability 
measures {P Ti ,...,t n }, where P Tll ..., TJV corresponds to the 
measure generated on O by the processes (Q , ■ ■ . , £t ) 
when the change in the TV-tuple process occurs at time point Ti, 
i = 1, . . . , TV. Notice that the measure Poo,...,oc corresponds 
to the measure generated on f2 by TV independent Brownian 
motions. 

Our objective is to find a stopping rule T that balances 
the trade-off between a small detection delay subject to a 
lower bound on the mean-time between false alarms and will 
ultimately detect min{ri, . . . , Tjv} H 

As a performance measure we consider 



(2) jif(T) = 
sup essup E Tl) 



{(T-ri A... At n })+\T TiA ... Atn } 



where the supremum over t% , . . . , tn is taken over the set 
in which min{ri, . . . , tn} < oo. That is, we consider the 
worst detection delay over all possible realizations of paths 
of the A^-tuple of stochastic processes (Q , ■ ■ ■ , Q ) up to 
min{ri, . . . , tat} and then consider the worst detection delay 
over all possible A^-tuples {ti, . . . ,tn} over a set in which 
at least one of them is forced to take a finite value. This is 
because T is a stopping rule meant to detect the minimum of 
the N change points and therefore if one of the N processes 
undergoes a regime change, any unit of time by which T 

2 Due to the symmetry of Brownian motion, without loss of generality, we 
can assume that fi > 0. 

3 In what follows we will use n A ... A tn to denote minjri , . . . , tjv}- 



delays in reacting, should be counted towards the detection 
delay. 

The criterion presented in (O results in the corresponding 
stochastic optimization problem of the form: 



(3) 



subject to E 



{T} > 7- 



We notice that the expectation in the above constraint is 
taken under the measure Poo,...,oo- This is the measure gener- 
ated on the space ft in the case that none of the N processes 
(^ 1 ',...,^ JV ') changes regime. Therefore, iOC {T} is 

the mean time between false alarms. 

In the case of the presence of only one stochastic process 
(say {£( 1 ' > }), the problem becomes one of detecting a one- 
sided change in a sequence of Brownian observations, or a 
vector of observations ■ ■ ■ jCt^) with the same change 
points, whose optimal solution was found in [3] and [15]. 
The optimal solution is the continuous time version of Page's 
CUSUM stopping rule, namely the first passage time of the 
process 



(4) y\ 



(i) 



sup loj 

0<Tl<t 



dp T 



dPr, 



,(!) 



where 



(5) u\ 
and 
(6) 



(i) 



/is, 



(1) 



-yt, 



,(!) 



inf u 

0<s<t 



(i) 



The CUSUM stopping rule is thus 

(7) T v = inf{t >0 ;2/t (1) >!/}, 

where v is chosen so that {T u } = ^-/(V) = 7, 
f{v) = e v — v — 1 (see for example [4]) and 



with 



(8) 



J$(T v ) = E {T v } = -sK-v) 



The fact that the worst detection delay is the same as that 
incurred in the case in which the change point is exactly is 
a consequence of the non-negativity of the CUSUM process, 
from which it follows that the worst detection delay occurs 
when the CUSUM process at the time of the change is at 
[4]. 

We remark here that if the N change points were the same 
then the problem ® is equivalent to observing only one 
stochastic process which is now A^-dimensional. Thus, in this 
case, the detection delay and mean time between false alarms 
are given by the formulas in the above paragraph. 

Returning to problem (O, it is easily seen that in seeking 
solutions to this problem, we can restrict our attention to 
stopping times that achieve the false alarm constraint with 
equality [8]. The optimality of the CUSUM stopping rule in 
the presence of only one observation process suggests that a 
CUSUM type of stopping rule might display similar optimality 
properties in the case of multiple observation processes. In 
particular, an intuitively appealing rule, when the detection of 



min{n, . . . , tat} is of interest, is 7), = A . . . A Tf, where 
Tl is the CUSUM stopping rule for the process {Q^;t > 0} 
for i = 1, . . . , N. That is, we use what is known as a multi- 
chart CUSUM stopping time [10], which can be written as 



(9) T h = 
where 

0) i 
Vt = SU P lo gj», 

0<Ti<t «"oO 



inf |i > 0; max-j^ 1 ^, 
dP T 



■Vi 



(AO 



}>h}, 



F, 



and the P Ti are the restrictions of the measure P Tl! ... lTN 
C[0,oo). 

It is easily seen that 



to 



J$\T h ) = E 



0,oo, . . .,oc 



00,0,00, . . .,oo 



,o ■ 



This is because the worst detection delay occurs when at 
least one of the N processes does not change regime. The 
reason for this lies in the fact that the CUSUM process is 
a monotone function of /i, resulting in a longer on average 
passage time if /j, = [9], That is, the worst detection delay 
will occur when none of the other processes changes regime 
and due to the non-negativity of the CUSUM process the worst 
detection delay will occur when the remaining one processes 
is exactly at 0. 

Notice that the threshold h is used for the multi-chart 
CUSUM stopping rule (O in order to distinguish it from v 
the threshold used for the one sided CUSUM stopping rule 
0. 

In what follows we will demonstrate asymptotic optimality 
of © as 7 — > oo. In view of the discussion in the previous 
paragraph, in order to assess the optimality properties of the 
multi-chart CUSUM rule (O we will thus need to begin by 
evaluating Bo,oo,...,oo {Th} and .,«, {T h }. 



Since the processes £. 



1, . . . , N, are independent it 



is possible to obtain a closed form expression through the 
formula 

(10) 

,oo,...,oo {T h } = / P Q ,oo,...,oo (Th > t) 

Jo 

T*0 , oo , .... oo ({Tl > *} n . . . n {if > t})dt 



/ PoiT^^lP^T^t^^dt. 
Jo 



Similarly, 



/•OO 

(11) E 00l ... t00 {T h } = / [PUTl >t)] N dt, 

Jo 

where {T^ > t} = {sup 0<s<t < h}. In other words, the 
evaluation of ( fTQb and ( fTTT is possible through the probability 



density function of the random variable sup 0<s<t ys for 
arbitrary fixed t which appears in [6], 

In order to demonstrate asymptotic optimality of (O we 
bound the detection delay of the unknown optimal 

stopping rule T* by 



(12) x {'//.} > J { m\T*), 

where h is chosen so that 

(13) /.\ = 7- 

It is also obvious that j¥?'(T*) is bounded from below by 
the detection delay of the one CUSUM when there is only one 
observation process, in view of the fact that 

supn,...,™ essupS n ,... jTJV | (T — 17 A ... A t n ) + \T TiA ... Atn } 

> su Pri essup E T1 { (T - n) + \T Tl } . 

The stopping time that minimizes 

sup T1 essup .E n |(T-Ti) + |^ 7 Tl | is the CUSUM stopping 
rule T v of (O, with v chosen so as to satisfy 

(14) E X {T V } = 7. 

We will demonstrate that the difference between the upper and 
lower bounds 

(15) h u x v {T h \ >jVp(T*) >E Q {T„}, 

is bounded by a constant as 7 — > 00, with h and v satisfying 
(fT3] > and (fl4] >. respectively. 
Lemma 1: We have 



> 



E { 
(16) 



0,00,. ..,00 i-^ftj o 



log 7 + log -l + o(l) 



as 7 — > 00, 

Proof: Please refer to the Appendix for a sketch of the proof. 
Moreover, it is easily seen from (JHJ that 



(17) E {T U } 



log7- 



Thus we have the following result. 



Theorem 1: The difference in detection delay of the 
unknown optimal stopping rule T* and the detection delay of 
Th of (0 with h satisfying (fT3l is bounded above by 

A (log N), 
M 

as 7 — » 00. 

Proof: The proof follows from Lemma Q] and ( TTTb . 

The upper and lower bounds on detection delay for the 
optimal stopping rule, when /.i is |, 1 and 2, in the case that 
N = 2 are shown in Figure [2] 

The consequence of Theorem Q] is the asymptotic optimality 
of (0 in the case in which all of the information becomes 
directly available through the filtration {J~t} at the fusion 
center. We notice however that this asymptotic optimality 
holds for any finite number of sensors N. 

We now discuss the implications of the above result in 
decentralized detection in the case of one shot schemes. 



The upper and lower bounds on the detection delay (DD) for the optimal stopping rule 



llupei Bound 
Lowef Bound 



?i:<j'.j 
T 





(b) fj, = 1 



(C) /i : 



(a) M = 5 

Fig. 2: (Left) Case of /i = |. (Middle) Case of /Z2 = 1. (Right) Case of /i = 2. (Note that the differences between upper and 
lowers bounds are all bounded as 7 increases.) 



III. Decentralized detection 

Let us now suppose that each of the observation processes 
{Q } become sequentially available at its corresponding sen- 
sor Si which then devises an asynchronous communication 
scheme to the central fusion center. In particular, sensor Si 
communicates to the central fusion center only when it wants 
to signal an alarm, which is elicited according to a CUSUM 
rule Tjp of 0. Once again the observations received at the N 
sensors are independent and can change dynamics at distinct 
unknown points tj. The fusion center, whose objective is to 
detect the first time when there is a change, devises a minimal 
strategy; that is, it declares that a change has occurred at the 
first instance when one of the sensors communicates an alarm. 
The implication of Theorem[T]is that in fact this strategy is the 
best that the fusion center can devise and that there is no loss 
in performance between the case in which the fusion center 
receives the raw data . . . £ t (A °} directly and the case in 

which the communication that takes place is limited to the 
one shown in Figure Q] To see this, the detection delay of the 
stopping rule T h = T^A. . .AT^ N] is equal to E 0<OOi ... tOO {T h } 
when Si is the one that first signals an alarm, -E 00 ,o,...,oo'P\} 
when 52 first signals and so on all of which are equal due 
to the assumed symmetry in the signal strength p received at 
each of the sensors 5, when a change occurs. The mean time 
between false alarms for the fusion center that devises the rule 



if A 



A T[ N) is thus E a 



,c {Th}. But Theorem 
Q] asserts that this rule, namely Th, is asymptotically optimal 
as the mean time between false alarms tends to 00 in the 
centralized case for any finite N. In other words, the CUSUM 
stopping rules T^, T^\ Tf are sufficient statistics (at 
least asymptotically) for the problem of quickest detection of 
©. 



thresholds {hi} should be chosen so that 



0,00,. ..,00 



{T c } — i?oo,0,c 



\Tc } — • • • — ■E o,...,oc,o{^C }) 



where Tq = Xj^ A ... A T^ n . 

A further interesting extension corresponds to the case of 
correlated sensors. To demonstrate this case let us begin by 
assuming that N = 2. This case corresponds to but with 

(18) £{w t (1) u4 2) j = pmhx{s,t} Vs,i>0, 

This case becomes significantly more difficult because of 
the presence of local time in the dynamics of the process 
max{y t , y t }. Nevertheless, it is possible to derive a for- 
mula for the expected delay of Th under the measure Poo, 00 ■ 
This expression is given by 

{Th} = \(e h -h-l) 



E 



2(l-p)£ j^V 



1)%. 



(i) 



yi 2) )ds\, 



(19) 



where S denotes the Dirac delta function and the final term 
in this expression corresponds to the collision local time of 
the processes and weighted by the factor (e Vt — 1). 
The difficulty in the use of expression ( fl9l l is the fact that as 
p changes, the expected value of the collision local time term, 
which is the last term in ( fl9l ), also changes. Moreover, the 
expression for the first moment of Th becomes significantly 
more complicated under the measure Po.oo- 



IV. Possible extensions 

An interesting extension corresponds to the case in which 
the signal strengths p are different in each sensor after the 
change. That is, after the change the signal in Si is fa with 
A*i 7^ A* 2 7^ • • ■ 7^ A*JV- in this case, it is not clear what 
the optimal choice of thresholds is, but it is possible that the 
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V. Appendix 

As an illustration for general case, let us prove the result 
for N = 2. 

We begin by deriving expressions for £0,00 {Th} and 
£-00,00 {Th} by using the results in [6]. For all h > 2, we 
have 



£0,00 {Th} = 

32 ^ sin 3 fa sin 3 @j 

7? 



sin fa sin 6j cos q 

cos fa)(6j — sin f?j cos 0j) cos 2 fa 



cos 2 fa cos 2 9j 



{Th} = 



e~ h sin 3 0i sin 3 1 



32 x ^ 

p? .4^ (0i — sin 0i cos 6i)(9j — sin 9j cos 6j ) cos 2 Oi + cos 2 0j 



G4 



1 sinh 77 



E 



sin 3 Oj cos 2 1 



cosh 77 



H 2 sinh 77 cosh 77 - 77 Oi - sin Oi cos t cos 2 Oi + cosh 77 

16 e"' 1 sinh 6 77 cosh 2 77 
fi 2 (sinh 77 cosh 77 — ?i) 2 

= S 4 (h) + S 5 {h)+Se(h), 



where 



tan <p n 
tan 6 n 
tanh 77 



< 



-9„ > 
h 

2 

^ >0 - 



The idea then is show Si(h), 82(h), 84(h) and 85(h) 
converge to zero, and examine how 82(h) and S5Q1) behave 
as h — * 00. In the following paragraphs we shall analyze 
these in the order S 6 (h) — > S 3 (h) — > S 2 (h) — * S 5 (h) — » 
S 1 (h)^S i (h). 

First notice that for large h, 77 is large and close to 4. 
Moreover, 



(20) 



2 2 "-' 1 = 1- 4^-^ + 0(6-^). 



This can help us to compare 77 with h. 
For S 6 (h), 

Se(h) = —er h sinh 4 77 (l — 7/ sinh 1 7/ cosh 1 77) 
fj, 

= \e- h [e 4r ' + (877-4) e 2 " + (e")l 

= -L [e^e 4 "- 2 ' 1 + (877 - 4)e 2 "-' 1 + o( e - r ')] , 
/_r 

by (l20l i the first term is 

e h e in-2h = e "( e 2 ^) 2 = e ' l [i_4 ??e - 2 '; + ( e -3^]- 
= e h [l - 87;e- 2r ' + o(e-^)} 
= e h - 8r]e h - 2n + o(e _ ") 
= e h - 877 + o(e _ "), 



- cos 2 



M 2 ^ (4>i - sin fa 

32 sinh 3 77 sm3 fa cos2 

/j 2 sinh 77 cosh 77 - 77 fa - sin fa cos fa cos 2 fa + cosh 2 77 



cos 2 fa 



J and the second term is 

(877 - A)e 2ri - h 



32 sinh 3 77 sm3 4>i cos2 4>i 

ft 2 sinh 77 cosh 77 — 77 fa — sin fa cos <j> 



(877-4) [l + o(e-")] 
877 - 4 + o(e- r '), 



so 



yi yi >^vjv5 yi 

= Si(h) + S 2 (h) + S 3 (h), 



(21) S 6 (h) 



e h -4 + o(e-2) 



, as h — > oo. 



For 83(h), also note that from ([H} and [6] we can write 

2 f°° 
—/(-ft) = / P Q {T h >t)dt 
M Jo 



(22) 



16 h 



re 2 



E 



sin 3 cos 2 0j 
f>i — sin cos < 



from which we obtain 

32 sinh 3 77 



S 3 (h) 



y? sinh 77 cosh 77 — 77 



E 

i>i 



sin 3 & cos 2 1 



sin 0j cos ( 



(23) 



2 



[l + o(e-")] (ft + e~' 1 - 1) 



/i - 1 + o(e~2 ) 



as ft — > 00. 



To bound 82(h) and S^ft) we need the following, 
Result 1: Suppose < p < 1. Then, for all positive solutions 

{cki}i>i t0 tne equation tana; = px (tana: = — px, resp.), we 

have 



(24) 



lim y 



I sin 3 ai I cos 2 Oj 1 



o^o+ ' — ' aj — sin a; cos a; 

z>l 



This suggests that, asymptotically, as h — > 00, 



sinh 77 



sinh 77 cosh 77 — 77 



< 

1 

- H 

7T 



i + (l) 

7T 



0(1) 



sin 4>i cos </>j cos 2 ^ + cosh 7? 

sinh 3 77 
cosh 2 ?7(sinh77COsh77 — 77) 

— — J^- ( 1 — 7? sinh -1 77 cosh -1 77) 
cosh 3 77 v 1 

= 0(e-i), 



from which we obtain 



(25) 



I $2 (ft) I = — jO(e 2 ), as ft — > 00. 



Similarly, 

e~ h sinh 3 77 

sinh 77 cosh 77 - 77 B l - sin 0j cos 0j cos 2 0, + cosh 2 77 



E 



sin 



cos 



cosh 77 



< 



1 



OK 



e h — — !— — (1 — rj sinh rj cosh 77) 1 



cosh 77 



0(e-»), 



so 
(26) 



|5 5 (ft) I 



r 0(e 2 ), as ft — > 00. 



Consequently, 

(28) |5i(ft)| 

Similarly, 



o(l), as ft — > 00. 



(29) 



1^4 (ft) I = — o(e h ), as ft -> 00. 



Finally, from ED, (E3, (E2]>, <ES> and (EH we obtain 



(30) E 0jOO (T h ) = -j [ft - 1 + o(l)] , as ft -> oo, 
and 

(31) £00,00(2*.) = -1 [e' 1 -4 + o(l)] , as ft ^00. 

And for ft and 7 satisfying ( TT3b , we have asymptotic results 
(O with N = 2. 

Now let us prove the two results we used in the above. 
Result 1: 

Proof: For any on £ ((•£ — i)7T, (« + |)tt) such that tana, = 
±pai, (0 < p < 1), we have 



I sin 3 cti I cos 2 a* 



p 3 a 2 



cti — sinajcosai (1 + p 2 af) 3 / 2 [(1 =F p) + P 2 ct 2 ] 
P P 



< 



(l+p 2 a 2 f/ 2 



< 



l+p 2 (i-\) 2 k 2 



3/2 ' 



Thus 

E 



sin 3 cti I cos 2 cti 



< 



E 



o.i — sin a, cos a.i 
i>i 1 11 i>x 



< 



pdx 



r 9 n3/2 

[l+p 2 (t-|) 7T 2 

cfot 



tt J* (1 + p 2 a; 2 ) 3 / 2 7r 7_ p i (1 + u 2 ) 



\3/2 



rZi/ 



= — , as p — > + . 



To handle the double sum in Si (ft) and Si(h), we need 
Result 2: Suppose < p < 1, {ai}i>i are all positive 
solutions to the equation tana; = px, and {/3j}j>i are all 
positive solutions to equation tana; = px (tana; = —px, 
resp.), then 

>— ^ sin 3 cti sin 3 f3j cos 2 an cos 2 /3j 

i j>l ^ ~~ S111 a * C ° S a *-"-^ ~~ S111 ^ C0S ^ ) C ° S2 ai C0S2 

(27) ->-0, asjj^0+. 



(1 + 7J 2 ) 3 / 2 7T 



Result 2: 

Proof: For simplicity let us denote the (i, j)-term in the sum 
by a>ij(p). As in the last proof, a little computation would 
give us 

\ a i,j(p)\ = Ip(pui,pf3j) -P 2 , 
where I p (x,y) (0 < p < 1) is the function 

I P (x,y) = 



v/(l + a; 2 )(l + y 2 )(2 + a; 2 + y 2 )(l + ^)(1 + ^) ' 

Clearly, 7 p (-, •) is (uniformly in p, < p < 1) bounded above 
by the L X (R 2 ) function B(-,-), which is defined as 

B(x,y) 



^/(l + x 2 )(l + y 2 )(2 + x 2 + y 2 ) 



We have two steps to finish our proof: 



(a) Urn \ 



") 2 



I (x 1 y)dxdy, 



(b) lim V 
p-o+ 4K 



Oi,j(p) ' 1 {a i ,j(?))>0} 



— / / I (x,y)dxdy. 



1 

Let us start from (a). Given any e > 0, we can find a 
constant M > such that, for Rm = {{x,y) : min(x, y) > 
M} and all < p < 1, 

{Ip is decreasing in both a; and y in Rm, 
^ I Ir m 7 f - / Jh m 2/)^ ^ I- 

Because of this, for any < p < 1, the "tail" sum 

£ 1 

min(pa.i ,pflj ) >M+p7r 

(32) 



i(p)\ < -% / / I P {x,y)dxdy < 



where we define cij,j(0) = for all (z, j). 

On the other hand, as p goes to zero, the function I p will 
converge uniformly in [0, M} 2 to Jo- So all the terms \a,ij(p)\ 
in the "head" sum are uniformly very close to Io {pcti , p(3j ) -p 2 , 
the sum of which, multiplied by tt 2 , is a Riemann sum of the 
function Io(x,y) over the region [0, M] 2 , and will converges 
to the Riemann integral of Io over [0, M} 2 as p turns to zero. 
In other words, for small p, there exists 



max(pa^ .pflj ) < M +p7r 



/(p)l 



l 



M /.M 



I (x 1 y)dxdy 



(33) 

By ( f32t and ( T33l l. we have 



lim y 



(P)l 



1 



;+)2 



I (x,y)dxdy 



< 



< e. 



(34) 



Now let e goes to zero we are done with (a). 

The proof of (b) is similar. Note that the signs of the 
a i,j(p)' s can be represented by (— or (— and in 
each rectangle [2(i — l)pw, 2ipTt] x [(J — l)pw,jpTr], > 1), 
either a2i-i,j (p) or a 2i,j (p) is positive. With the same constant 
M chosen as above, for the sum of all positive aij(p)'s such 
that max(pai,p(3j) < M+pir, we can use the same argument 
as before, to show that for small p, 



Thus (b) is proven because both the tail integral and the tail 
sum are negligible due the way to choose M. ■ 
In the N CUSUMs case with N > 2, the calculation 
is similar: both of the main terms in £7o,oo,...,oo {Th} and 
in £Joo,oo,...,oo {Th} are the terms with highest degree in 
[sinh 3 r)/ (sinh 77 cosh r\ — r))] . With ( f20b we can get they are 

%[h-l + o{e-%)], 



(36) 
and 
(37) 



e h + {N- 2)h + (2 - 37V) + o(e - T) 



N/i 2 
respectively. 

We can prove that all other terms converge to zero as h 
goes to infinity. 

With a generalization of Result 1 to n dimensional trigono- 
metric sums and integrals for all n > 1, we are able to 
deal with most terms in the expansion of the expectations, 
because those bounded trigonometric sums are multiplied by 
expressions of negative exponential order in h. 

There is only one term (in £?o,oo,...,oc {Th}) which cannot 
be proven to converge to zero in this manner. We need to 
prove the sum involved there, which is ( f38l ) at the top of 
the following page, converges to zero as h goes to infinity. 
We can follow the proof of Result 2 to get the result. To be 
more precise, denote p = |, and the term in above sum by 

N ?(p)\ < \af]{p)\, that can help 



(JV) 

us to control the "tail" sum 



a,-,- (p), then obviously. 



(40) 



min(p0i ,p0j ) > M+J37T 



4?' 



<«-.• (p)i, 



where M is chosen as in the proof of Result 2. On the other 
hand, 



(41) 



a\ N Mp)\=li N Hp^ P J )-p 2 



where Ip N ^ is the function defined in d39l ). The function Ip 
uniformly converges to Io as p goes to infinity in the domain 
[0, M} 2 , since pr\ — > 1 as p — > + . As a result, the "head" 
sum converges to the same double integral as the one in j33l 
or (|35l l, so we are done! 

Finally, by (|36l l and (|37l i, we can derive asymptotic formula 
([Tol l with h and 7 satisfying (fT3l . 



(AO 



' 2 "~ i a ijW ' I K, 3 (p)>o} 



(35) 



M 



Jo 



I (x,y)dxdy. 



(38) ]T 



sin 3 (pi sin 3 9j cos 2 <pi cos 2 9j 



j >1 {<t>i — sin fa cos fa) (9 j — s'm9j cos 9j)[(N — 2)(1 — 4^) cos 2 ^ cos 2 9j + cos 2 ^ + cos 2 9j] 

(39) l( N Hx v) = ^ - 

p y/(l + x 2 )(l + y 2 )[(7V - 2)(l-pV) + 2 + .x 2 + y 2 ](l + + ±#) 



